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A model-independent partial-wave analysis of the S'-wave component of the Kir system from 
decays of D + mesons to the three-body K~-k + -k + final state is described. Data come from the 
Fermilab E791 experiment. Amplitude measurements are made independently for ranges of K~tt + 
invariant mass, and results are obtained below 825 MeV/c 2 , where previous measurements exist only 
in two mass bins. This method of parametrizing a three-body decay amplitude represents a new 
approach to analysing such decays. Though no model is required for the S-wave, a parametrization 
of the relatively well-known reference P- and D-waves, optimized to describe the data used, is 
required. In this paper, a Breit-Wigner model is adopted to describe the resonances in these waves. 
The observed phase variation for the S-, P- and D- waves do not match existing measurements of 
I = i K~tv + scattering in the invariant mass range in which scattering is predominantly elastic. 
If the data are mostly I = this observation indicates that the Watson theorem, which requires 
these phases to have the same dependence on invariant mass, may not apply to these decays without 
allowing for some interaction with the other pion. The production rate of K~ir + from these decays, 
if assumed to be predominantly I = ^, is also found to have a significant dependence on invariant 
mass in the region above 1.25 GeV/c 2 . These measurements can provide a relatively model-free 
basis for future attempts to determine which strange scalar amplitudes contribute to the decays. 
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I. INTRODUCTION 

Kinematics and angular momentum conservation in decays of ground state, heavy-quark mesons to three pseu- 
doscalars strongly favor production of S-wave systems. These decays have therefore been regarded as a source of 
information on the composition of the scalar meson (spin-parity J p = + ) spectrum. Extracting this information 
has, however, been done in model-dependent ways that can influence the outcome. For the di-mcson subsystems, 
vector and tensor resonances are relatively well-understood, but, as larger samples of D and B meson decays become 
available, the correct modelling of the S-wave contributions becomes an increasingly important factor in the task of 
obtaining satisfactory fits to the data. 

Analyses typically use an isobar model formulation in which the decays are described by a coherent sum of a 
non-resonant three-body amplitude NR, usually taken to be constant in magnitude and phase over the entire Dalitz 
plot, and a number of quasi two-body (resonance + bachelor) amplitudes where the bachelor particle is one of the 
three final state products, and the resonance decays to the remaining pair. It is assumed that all resonant and NR 
processes taking part in the decay are described by amplitudes that interfere, and have relative phases and magnitudes 
determined by the decay of the parent meson. In cases where all three decay products are pseudoscalar (P) particles, 
angular momentum conservation requires that the resonances produced are scalar (5- wave), vector (P-wave), etc . 
For D mesons, decays beyond D-wave are highly suppressed by the angular momentum barrier factor and can be 
neglected. 

Within this formalism, the decays D + — > 7r~7r + 7r + and D + — > K~it + ii + pj were once thought to require very large, 
constant NR amplitudes 0- [3- [1 Using larger samples, the Fermilab E791 collaboration found that a satisfactory 
description of these decays requires more structure. By including S'-wave isobars, <r(500) — > 7r + 7r _ in 7r~7r + 7r + [j| and 
k(800) — > K~n + in K~it + it + |(|, a much-improved modelling of the Dalitz plots was achieved, and the need for a 
constant NR term was much reduced in each case. 

The FOCUS collaboration, using an even larger sample of D + — > 7r~7r + 7r + decays, found an acceptable fit 
using a ifmatrix description of the S'-wave with no cr(500) pole. However, a parametrization of the NR background 
was required to achieve an acceptable fit. The BaBar and Belle collaborations H, 0, E|j with the measurement of 
CP violation parameters in B~ — ► D a (^ K°ir + ir~)K~ decays as their primary goal, introduce ct(500) and another 
ct(1000) isobar in order to obtain an acceptable description of the complex amplitude for the D° Dalitz plot. 

The important issue of whether scalar particles cr(500) and k(800) exist is not convincingly settled. Further 
observations of these isobars were recently reported in tt~tt + and K~ir + systems from J/ij) decays |Tllll2j. However, 
these results were modelled on variations of the simple Breit-Wigner form for the states adopted in the cases cited. 
Quite different descriptions are probably required, since such forms are seen to contain poles below threshold fili |. In 
a recent publication [TJ| the E791 data on D + — ► 7r~7r + 7r + and on the D + — ► K~n + Tr + decays discussed here, are 
re-fitted using input from calculations of irir and Kit scattering that include constraints of chiral perturbation theory, 
and that find both a and k poles 0,^1- The fits obtained yield similar \ 2 P er degree of freedom to those in Rcfs. 
[j| and where Breit-Wigncrs were used, but each resonance is considerably wider. 

Ultimately, a less model-dependent analysis of the data should help resolve the issue of the a and the k. 

Model-independent measurements of the energy dependence of these S-wave amplitudes, particularly in the low 
invariant mass regions, where confusion is greatest, is therefore an important experimental goal. Such a Model- 
Independent Partial Wave Analysis (MIPWA) is reported here for the K~ir + system produced in D + — > K~tt + it + 
decays. One earlier measurement has been made for tt~tt + systems from D + — ► tt~7t + tt + decays |l7j . in which the 
"amplitude difference" (AD) method was employed. This method can only be used when there exists a region 
of the Dalitz plot that can be described by the sum of a single resonance and an S-wave amplitude that is to be 
measured. Interference of the resonance with this S-wave introduces an asymmetry in the distribution of the other 
invariant mass combinations that can be measured at different values of invariant masses in the band. As there is no 
such region in the Dalitz plot for the data reported here, this method is not used. 

For the K~tt + system, the best results of an MIPWA currently available come from the LASS experiment [l^, in 
which K~ir + scattering was studied for invariant masses above 825 MeV/c 2 . Below 825 MeV/c 2 , measurements have 
been made for the mass bins 770-790 McV/c 2 [2C| and 700-760 McV/c 2 [2ll|. though with less precision. Information 
on the Kir S-wave amplitude near or slightly above the JC*(892) has been extracted by the BaBar collaboration in 
studies of B decays to J/ipKir [23] , and by FOCUS in semi-leptonic D decays to Kirlv [2^ , the low mass region has 
not been covered in either case. 

In this paper, we describe an MIPWA in the mass range from K~n + threshold up to 1.72 GeV/c 2 , the kinematic 
limit for decays of D + mesons to K~tt + tt + final states. The amplitudes obtained for the S-wave require no assumptions 
about its dependence on invariant mass, though they do rely on a model for the relatively well-understood P- and 
D-waves. As such, they should provide an unbiased input for comparisons with theoretical models for scalar states. 

This paper is organized as follows. In the following section we present the data sample. Next we describe the 
method used to extract complex amplitudes from the S-wave K~ir + system in a way that does not require a model 
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for its dependence on invariant mass. In Section llVl this is applied to the sample of D + — > K~ir + ir + decays. The 
amplitudes obtained are then compared, in Section [V] with the S-wave amplitude derived from the Breit-Wigner 
isobar model fit that best represents the data. This model, applied to these data, was presented in Ref. @, and 
includes a k(800) isobar. In Section IVll the results of the MIPWA are compared with amplitudes measured in K~ir + 
elastic scattering, and with the expectations of the Watson theorem |24j , whose applicability to weak hadronic decays 
has not previously been tested. Systematic uncertainties are discussed in Sec. I VIII Finally, some conclusions are 
drawn. In Appendix^we point out limitations, ambiguities and other technicalities inherent in this kind of analysis. 



II. THE E791 DATA 



The analysis is based on a sample of D + — > K~ir + ir + candidates from Fcrmilab experiment E791. The experiment 
is described in detail in Ref. [2^|. The same sample is used in this paper as the one described in Ref. [(|, where an 
isobar model fit to these data was described. The selection process used in obtaining the sample is outlined below, 
but more details are given in Ref. @ . 

In this paper, s is used to denote the K~ir + squared invariant mass. Where it is important to distinguish, the 
two pions (and their corresponding s values) arc labelled, respectively, t\\ and 7rJ (s A and s B ). A clear peak in the 
K~ir + ir + invariant mass M distribution is observed with 15,079 events in the mass range 1.810 < M < 1.890 GeV/c 2 , 
of which 94.4% arc determined to be signal. The major sources of background are incorrect three-body combinations 

(3.58%), and reflections of Df — > (f>TT + and Df — > K K + decays (1.75% and 2.61%, respectively) in which a K + is 
incorrectly identified as a ir + . The probability density function (PDF) for these backgrounds over the Dalitz plot is 
obtained from events in the sideband region of the K~tt + tt + invariant mass distribution and, for the second and third 
sources, from a large sample of Monte Carlo (MC) simulated events. An appropriately weighted combination of these 
three backgrounds is determined from their distributions in M . The efficiency for reconstructing the D + decays (the 
signal) is also determined from the MC events. It is described in this paper by a function e(s A , s B ). 



A. E791 Dalitz Plot 



The symmetrized Dalitz plot for this sample is shown in Fig. n wncrc s a is plotted vs. s B (and the converse). A 
horizontal (and the symmetrized vertical) band corresponding to the presence of the P-wave AT* (892) resonance is 
clear. Complex patterns of both constructive and destructive interference near 1400 McV/c 2 due to either o-wave 
(1430), P-wave A*(1410) or D-wave AT|(1430) are also observed. A further contribution from the P-wave A* (1680) 
state is also present, as determined by fitting. This is difficult to see due to smearing of the Dalitz boundary resulting 
from the finite resolution in the three-body mass. 

All these resonances are well-established and are known to have significant K~ir + partial widths. Interference 
between resonances is evident in the regions of overlap. 



B. Asymmetry in the K n System 



One of the most striking features of the Dalitz plot is the asymmetry in each K*(892) band. In any given K~ir + 
mass slice, a greater density of events exists at one end of that slice than at the other. This asymmetry is also evident 
in the region closest to the if* (892) peak itself. This is most readily explained by interference with a K~ir + S'-wave 
component and clearly shows that these data can be used to infer the structure of the S'-wave amplitude, provided 
an adequate modelling of the remainder of the plot is possible. 

This asymmetry, a, depends on the distribution of the hclicity angle, 9, the angle between K~ and 7r+ in the A' _ 7r+ 
rest frame. It is defined |2g as 

N CO s8>0 — A cos 0<o 

a= N li W 

where N is the efficiency-corrected number of events in the indicated regions of cos 9. In Fig. [3 a is plotted as 
a function of the AT*(892) Breit-Wigner (BW) phase <fi Blv — tan -1 [moT/(m,Q — s)], where the peak mass mo = 
896.1 MeV/c 2 and the mass dependent width T = 50.7 MeV/c 2 at the peak mass. A change in the sign of a occurs 
when (f) BW ~ 56°, at an invariant mass below the A'*(892) peak. We note here that, in Ktt elastic scattering 0], a 
is observed to reach zero at <fi BW w 135.5°, a mass above the K*(892) peak. Evidently there is a ~ —79° shift in s-p 
relative phase in this D + decay relative to that observed in Ktt elastic scattering. 
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FIG. 1: Dalitz plot for D — > K~tt a -k b decays. The squared invariant mass s B of one K~tt combination is plotted against 
s A , the squared invariant mass of the other combination. The plot is symmetrized, each event appearing twice. Lines in both 
directions indicate values equally spaced in squared effective mass at each of which the 5*-wave amplitude is determined by the 
MIPWA described in section ITTT1 Kinematic boundaries for the Dalitz plot are drawn for three-body mass values M = 1.810 
and M = 1.890 GeV/c 2 , between which data are selected for the fits. 



III. FORMALISM 



A. K n Partial Wave Expansion 

The Dalitz plot in Fig.^is described by a complex amplitude Bose-symmetrized with respect to the identical pions 
7rj~ and ~k% : 

A = A(s A ,s B ) + A(s B ,s A ). (2) 

Considering the simplest, tree-level quark diagrams, iso-spin 1 = 1/2 K~n + systems are most likely to be produced. 
The contribution of the t7 + tt + amplitude to these decays is not expected to be significant, coming mostly from re- 
scattering processes. To test this, data are taken from measurements of n + p — > TT + n + n reactions |27| in which the 
phase of the t7 + tt + amplitude was found to vary slowly, assuming it to be elastic, from zero at threshold to about 
—30° at 1.45 GeV/c 2 , the upper range of the measurements. No evidence for isospin 1 = 2 resonances exists in this 
range. This amplitude is added to those in model C in Ref. It is found that the n + ir + contribution is, indeed, 
insignificantly small (0.7 ± 0.4)%. 

The amplitude A is therefore written as the sum of K~-k + partial waves labelled by angular momentum quantum 
number L, 

A{s A7 s B ) = J2(-2pq) L P L {cos9)x 

L=0 

T L D (q,r D ) x C L (s A ), (3) 

corresponding to production of K~tt + systems with spin J = L and parity (— 1) L in these D + decays. In this analysis, 
the sum is truncated at L max = 2 since the D-wavc A'|(1430), as measured in reference 0, contributes only about 
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FIG. 2: The asymmetry a plotted vs. BW phase <f>B\v for the K*(892). These quantities are described in the text, a becomes 
zero at 4> B w ~ 56 degrees. 



0.5% to the decays. This is already small and higher partial-waves are expected to be even further suppressed by 
the angular momentum barrier. With no way to distinguish I = \ and / = | components in the K~ir + systems 
produced, their sum is measured in this paper. 

In Eq. @, p and q are momenta for the K~ and bachelor 7r+ respectively, in the K~~k\ rest frame. The cosine of 
the hclicity angle is then given in terms of the masses m K - (m 7r +) and energies E K - (E v +) of the K~ (7rj) in the 
A'~7r+ rest frame by: 



cos ( 



p-q 



E K -E n+ -(s B -m 2 K .-ml + )/2 



pq 



(4) 



This is the argument of the Legendre polynomial functions Pl- J 7 ^, is a form factor for the parent D meson which 
depends on q, L and on the D's effective radius r = r D : 



TO = ~(rq) 2 /12 

J* = [l + (rg) 2 ] _ * 

Tl = [9 + 3(rg) 2 + (r 9 ) 4 ]-^ 



scalar 
vector 
tensor 



(5) 



For L > 0, these form- factors are derived for non-relativistic potential scattering J28j. For L = 0, the Gaussian form 
in Eq. (0, suggested by Tornqvist j2|j to be a preferred way to describe scalar systems, is used. This form was used 
also in Ref. 

The Cl(s a ) are complex functions, and are the invariant-mass-dependent parts of the respective partial waves. 
They do not depend on the other Dalitz plot variable s B and are referred to in this paper as the K~tt + amplitudes. 
Provided that interactions between the K~-k\ system and the bachelor 7r+ can be neglected, the Cl(s a ) are related 
to the corresponding amplitudes, Tl(s) measured in K~ir + scattering experiments, by 



C L {s) = |C L (s)|e^W = 



r sV L {s)T L {s) 



P 



(6) 



where Vl(s), unknown functions, describe the K~n + production in each wave in the D decay process pfflj . These 
replace the K~Tt + coupling present in elastic scattering (proportional to the 2-body phase-space factor y/s/p and 
barrier factor p L ). 

The principal goal of this analysis is to measure Co(s), using all higher L contributions to the Dalitz plot as an 
"interferometer". This requires a model for Ci(s) and C^s), the reference P- and D- waves. 
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B. The Reference Waves 



As in previous analyses, a Breit-Wigner isobar model is used to describe the P- and D- waves. Linear combinations 
of resonant propagators W H , one for each of the established resonances R having the appropriate spin, and each with a 
complex coupling coefficient with respect to K*(892), B R = b R e l ^ R , are constructed. Three possible K~tt + resonances 
are included in the P-wave, but only one in the -D-wave in the invariant mass range available to these decays: 

Ci(s) = [Wif» (892) (s) + B K , (I4I0) Wir» (1410) (s) + 



B„ 



) W K «(i680)(s)] x F R (p,r R ), 



C 2 (s) = S K »(i4 30) >V K »(i4 3 o)(s) xT R (p,r R ). 



(7) 

(8) 



where T R is a form factor for the resonances in the iC _ 7r + system, required to ensure that the resonant amplitudes 
vanish for invariant masses far above the pole masses. It is assumed to have the same dependence on ccnter-of-mass 
momentum and angular momentum as the D form factor !F^, but to depend on a different effective radius r = r R . The 
coefficients in Eq. (J7J) have their origin in the K~ir + production process arising from D + decays, and are therefore 
treated as unknown parameters in the fits. 

Each propagator is assumed to have a Breit-Wigner form defined as: 



W„(a) = 



1 



m R — s — im R T(r R , s) ' 



where m R and F R are the resonance mass and width, and 

r(r H ,s) = 

where p R is the value of p when s = m R . 



( m H \ 






\V~sJ 




T^{p R ,r R )_ 



(9) 



(10) 



C. Parametrization of the 5-wave 



The goal is to define the 5-wave amplitude making no assumptions about cither its scalar meson composition, nor 
of the form of any 5*-wave NR terms. To this end, two real parameters are introduced 

ca; = |Co(sfc)| ; 7fc = <Po(s k ) (11) 

to define the amplitude Co(sfc) = Cke llk at each of a set of invariant mass squared values s = Sk (k = 1, N s ). A second 
order spline interpolation is used to define the amplitude between these points (sk, Cke %lk ) The ck and 7fc values 
are treated as model-independent parameters, and are determined by a fit to the data. 

To obtain the results in this paper, N s = 40 equally spaced values of are chosen. These are indicated by the 
lines drawn on the Dalitz plot in Fig. ^ Other sets of values for Sk are also used to check the stability of the results 
obtained. 



D. Maximum Likelihood Fit 



In this analysis, the 3-body mass M is not constrained to be that of the D + meson. The fits are therefore made in 
three dimensions (M, s A , s B ). A normalized, log-likelihood function is defined as 



C 



events 




(12) 



where P s and P b * are the normalized signal and background PDF's, respectively. 

Three backgrounds (i = 1, 2, 3), described in Sec.|nl are included incoherently in Eq. (|12f) . Each is considered to 
constitute a fraction /j of the event sample in the selected range 1.850 < M < 1.890 GeV/c 2 , and to be described by 
the PDF: 



P 1 



i{M)6i{s A ,s B ) 



(13) 
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This expression has a three-body mass profile Qi(M) and a distribution 9i{s A , s B ), with normalization rij, over the 
Dalitz plot. For the combinatorial background, the PDF is determined by events in a band of M values above the 
D + peak, while for the D s reflections it is determined from the simulated MC samples. 
The signal PDF is 

p = Q Q {M)e{ SAl s B )\A{ SA ,s B )\ 2 

Jds A ds B dM F(M)e(s A ,s B )\A(s A ,8 B )\2' [ ' 

in which Qq(M) describes the shape of the signal component in the K~n + ir + invariant mass spectrum, parametrized 
as the sum of two Gaussian functions, and e(s A , s B ) is the efficiency for reconstructing these events. The normalization 
integral extends over the entire Dalitz plot for each M in the selected range. 



E. Decay Channels and Branching Fractions 



The amplitude A(s A , s B ) in Eq. J3J can be written as a sum over the N ch possible decay channels of the D + : 

A(s A ,s B ) = ^2A k , (15) 
fc=i 

where Ak is the complex amplitude for the k th decay mode for decay to the K~tt + tt + final state through either a 
resonance, or through the whole set of possible 5-wave and NR states. The fraction, F k , is computed for each such 
mode: 

f DP \A k \ 2 ds A ds B 

" L P \EiM 2 ds A ds B - { > 

This is the definition most often used in the literature on three body decays. It guarantees that each F k is positive. 
Due to interference, however, the F k do not necessarily sum to unity. 



F. Parameters, Phases and Constants 



The log-likelihood, Eq. |Q , is defined by many parameters. By choice, a number of these are held constant in 
the fits. Parameters for the background models P£ and their fractions /* are determined by studies of data and of 
MC samples and are fixed. Masses and widths for well-established P- and D-wave resonances are also held constant 
at values listed in Table [I] These come mostly from the Review of Particle Properties (RPP) publication For 
the if* (1680) values appropriate for the state with known coupling to Kir observed in Kir scattering in the LASS 
experiment |l9l | are used. The form factor radii are fixed at r D = 5.0 GcV -1 and r H = 1.6 GcV -1 , values determined 
in Ref. to be those providing the best isobar model description for these data. Isobar coefficients B R and partial 
wave amplitude parameters Cj and ji are generally allowed to vary. 



TABLE I: resonance mass m R and width F R values used in the fits described in this paper. With the exception of if* (1680), 

parameters are as quoted in Ref. 13^1 . 

Resonance m R (MeV/c 2 ) F R (MeV/c 2 ) 

if* (892) 896.1 50.7 

if t (1410) 1414.0 232 

if t (1680) 1677.0 205 

if 2 * (1430) 1432.4 109 



Phases are defined relative to the if* (892) resonance. In all fits described here, the coefficient for the if* (892) is 
taken to be real and of magnitude unity, as explicit in Eq. 

Two sources of uncertainty in this method result from the parametrization of the P-wave, and from the fact that 
several local minima in the likelihood function exist. These limitations are discussed in more detail in Appendix lAl 
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IV. MIPWA OF THE A' tt + 5-wave 



The technique described in Section ITTT1 is applied to the data shown in the Dalitz plot in Fig.^ The P- and P-wave 
amplitudes defined as in Eqs. Q and (JHJ) arc chosen as reference waves. The 40 equally spaced values Sk, indicated by 
lines in the figure, are chosen. The 5*-wave magnitude and phase , and 7^., at each Sk, and the P-wave and P-wave 
couplings Bi are all determined by the fit. With all established vector and tensor resonances with masses and widths 
shown in Table there are 86 free parameters. 

It is confirmed that the contribution from A"* (1410) is negligible, as reported in Ref. and this is dropped from 
further consideration. The fit is made with the remaining 84 free parameters. The complex coefficients B t and the 
fractions F t for each of the resonances i included in the P- and P-waves are summarized in Table [n] |33j 



TABLE II: Fractions, magnitudes and phases for resonant and 5*-wave components from four fits to decays of D + mesons to 
K~ir + ir + described in the text. Fit labels are "MIPWA" for the fit, described in Sec. II VI where magnitudes and phases for 40 
A _ 7T + mass slices described in the text are free to vary. Systematic errors are included for this fit. "Isobar" refers to the fit 
described in Sec.[V] The fit labelled as "Elastic" is described in Section IVT1 



Channel 


Fit 


Fraction 


Amplitude 


Phase 






F% 


b 


(3 (degrees) 


A*(892)7r+ 


MIPWA 


11.9±0.2± 2.0 


1.00 (fixed) 


0.0 (fixed) 




Isobar 


12.6± 1.6 


1.00 (fixed) 


0.0 (fixed) 




Elastic 


12.8± 2.0 


1.00 (fixed) 


0.0 (fixed) 


Ai(1680)7r+ 


MIPWA 


1.2±0.6± 1.2 


1.63±0.4± 0.2 


42.8±16.3± 4.5 




Isobar 


2.1± 0.4 


2.18±0.2 


28. 2± 7.2 




Elastic 


5.0± 0.8 


3.15±0.3 


17.1± 7.5 


A 2 *(1430)tt + 


MIPWA 


0.2±0.1± 0.1 


4.31±1.0± 1.1 


-12.2±23.7± 16.8 




Isobar 


0.5± 0.1 


6.50±0.7 


-54.0± 7.4 




Elastic 


0.3± 0.1 


4.59±0.0 


-46.9±12.3 


Total S-wave: 












MIPWA 


78.6±1.4± 1.8 


EIPWA 


EIPWA 




Elastic 


79. 2± 1.1 






S-wave components: 










NR 


Isobar 


16.1± 5.3 


0.60±0.1 


-3.5± 9.1 


Kn + 


Isobar 


45.6±10.7 


1.71±0.2 


181. 3± 8.1 


A5(1430)tt+ 


Isobar 


12. 2± 1.3 


0.52±0.1 


47.0± 5.6 



The S'-wave phases jk (= </>o(sfc)) and magnitudes Cfc (= |Co(sfc)|) resulting from the fit are plotted, with error 
bars, in Figs.|3Ja) and (b), respectively. A significant phase variation is observed over the full range of invariant mass, 
with the strongest variation near the Aq(1430) resonance. The magnitude is largest just above threshold, peaking at 
about 0.725 GeV/c 2 , above which it falls. A shoulder is seen at the mass of the Aq(1430), after which the magnitude 
falls sharply to its minimum value just above 1.5 GeV/c 2 . 

The S-wave magnitudes Cfc obtained depend on the form used for T° D in Eq. ©. The products c^-P^, and phases ~fk 
are, however, independent of T a D . To simplify future comparisons, values for Cfc, T° D and 7^ for each invariant mass Sk 
are listed, with their uncertainties, in Table ITTT1 In the present analysis, the Gaussian form [2|j in Eq. JSJ for L = 
has been chosen. The values used for T° D at each sj, are also listed in Table ITTT1 

The magnitudes |Cl(s)| and phases 4>l{s) for the P- and P-wave amplitudes Cl{s) (L = 1, 2) are computed from 
Eqs. (|7J) and JSJ, using parameters for this fit from Tables HI and HTI Uncertainties in these quantities are also computed, 
using the full error matrix from the fit. Values, at each s, plus or minus one standard deviation are then plotted as 
solid curves, with shading between them, in Fig. [3J The P-wave phase and magnitude are shown, respectively, in 
Figs.EUc) and (d), and those for the D-w&ve in (e) and (f). 

To compare the fit with the data, MC simulated samples of events are produced in the three-dimensional space 
in which the fits are made. Events are generated with the distribution predicted from the signal and background 
PDF's defined in Eq. (|12fl . Parameter values from Tables ITT1 and IIIl| and the measured event reconstruction efficiency 
e(s A ,s B ), arc used in the simulation. These events are projected onto the two-dimensional Dalitz plot, and its one- 
dimensional invariant mass plots. Data are then overlayed for comparison. These plots are shown in Fig. 0] As a 
further comparison, the distributions of the hclicity angle 9 in the A' _ 7r + systems predicted by the fit are compared 
to the data. Fig. [S] shows moments for this angle, (dN / dm) (Pl (cos 8)}, uncorrected for acceptance. 

Qualitatively, agreement between the fit and the data is very good. Quantitative comparison is made using the 
observed distribution of events on the Dalitz plot. The plot is divided into rectangular bins. For each of these, the 
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FIG. 3: (a) Phases jk = <^o(sfc) and (b) magnitudes Ck = |Co(sfc)l of S-wave amplitudes for K~n + systems from D + — > K~-k + tv + 
decays with the amplitude and phase of the _ftT*(892) as reference. Solid circles, with error bars show the values obtained from 
the MIPWA fit described in the text. The effect of adding systematic uncertainties in quadrature is indicated by extensions 
on the error bars. The P-wave and D-wave phases are plotted in (c) and (e) and their magnitudes in (d) and (f), respectively. 
These curves are derived from Eqs. J2J and JHJ, respectively, evaluated with the parameters and error matrix resulting from 
the MIPWA. Curves appear as shaded areas bounded by solid line curves representing one standard deviation limits for these 
quantities. In all plots, the dashed curves show one standard deviation limits for the predictions of the isobar model fit described 
in Sec. [V] These curves are computed in the same way, using Eq. 11711 in addition to Q and JSJ with parameters and error 
matrix from the isobar model fit. 



normalized residual, (n/ — n )/cr, where n a is the number of events observed, nj is the number predicted by the fit 
and a = ^Jn / + n a is the uncertainty in rif — n , is computed. The expected population in each bin, n/, is computed 
by numerical integration of the PDF in Eq. 1121 Neighboring bins are combined, where necessary, to ensure that 
n/ > 10. The normalized residuals, plotted as an inset in Fig.^d), are combined to obtain x 2 /NDF where NDF is the 
number of bins less the number of free parameters in the fit. These values, and the probability for obtaining them, 
are tabulated with the optimum log-likelihood value from the fit, given in Table llVl 



V. COMPARISON WITH AN ISOBAR MODEL FIT 



It is interesting to compare the results from the MIPWA with those reported in Ref. 6] which came from a Breit- 
Wigncr isobar model fit. In this fit, the S'-wave was modelled as a sum of isobars with Breit-Wigner propagators for 
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TABLE III: Magnitude c and phase 7 of the K~tt + S-wave amplitude determined, at equally spaced masses, by the MIPWA 
fit described in the text. The magnitudes assume a real form-factor F" for the D + meson. Values for this form-factor are given 
for each mass value in the table. Two further mass values, used in the fit as a result of the finite resolution in 3-body invariant 
mass, are not included in the table. They lie, respectively, at and above the kinematic limit for D + decay to this final state. 
Statistical and systematic uncertainties are assigned to each magnitude and phase. 



(GeV/c 2 ) 








c 

(GeV/c 2 )- 2 










7 

(degrees) 






0.672 


0.26 


8.37 




0.73 


± 


0.62 


-102 


± 


5 


± 


3 


0.719 


0.27 


9.04 


± 


0.59 


± 


0.89 


—96 


± 


5 


± 


3 


0.764 


0.29 


7.82 


± 


0.54 


± 


0.89 


—73 


± 


g 


± 


4 


0.807 


0.31 


7.42 


± 


0.43 


± 


0.57 


—77 


± 


7 


± 


5 


0.847 


0.33 


6.47 


± 


0.30 


± 


0.46 


—60 


± 


5 


± 


6 


0.885 


0.34 


5.57 


± 


0.31 


± 


0.07 


—54 


± 


6 


± 


5 


0.922 


0.36 


5.90 


± 


0.46 


± 


0.09 


—68 


± 


8 


± 


7 


0.958 


0.38 


6.17 


± 


0.52 


± 


0.01 


—72 


± 


10 


± 


9 


0.992 


0.40 


4.87 


± 


0.35 


± 


0.19 


—41 


± 


12 


± 


10 


1.025 


0.42 


4.42 


± 


0.28 


± 


0.09 


—43 


± 


11 


± 


5 


1.057 


0.44 


4.02 


± 


0.26 


± 


0.01 


—38 


± 


12 


± 


5 


1.088 


0.46 


3.74 


± 


0.19 


± 


0.11 


—22 


± 


10 


± 


4 


1.118 


0.49 


3.81 


± 


0.19 


± 


0.13 


—29 


± 


9 


± 


4 


1.147 


0.51 


3.16 


± 


0.14 


± 


0.13 


_3 


± 


g 


± 


4 


1.176 


0.53 


3.21 


± 


0.15 


± 


0.13 


_11 


± 


7 


± 


3 


1.204 


0.55 


2.86 


± 


0.14 


± 


0.32 


_3 


± 


7 


± 


3 


1.231 


0.58 


3.11 


± 


0.15 


± 


0.13 


-3 


± 


6 


± 


2 


1.258 


0.60 


2.92 


± 


0.15 


± 


0.24 


8 


± 


6 


± 


3 


1.284 


0.62 


2.80 


± 


0.16 


± 


0.18 


11 


± 


6 


± 


2 


1.310 


0.65 


2.77 


± 


0.17 


± 


0.12 


11 


± 


5 


± 


2 


1.335 


0.67 


2.83 


± 


0.17 


± 


0.20 


22 


± 


5 


± 


2 


1.360 


0.69 


2.73 


± 


0.19 


± 


0.31 


31 


± 


4 


± 


2 


1.384 


0.71 


2.29 


± 


0.20 


± 


0.25 


30 


± 


5 


± 


2 


1.408 


0.74 


2.38 


± 


0.23 


± 


0.01 


46 


± 


4 


± 


2 


1.431 


0.76 


2.05 


± 


0.28 


± 


0.08 


55 


± 


4 


± 


2 


1.454 


0.78 


1.59 


± 


0.25 


± 


0.07 


64 


± 


6 


± 


4 


1.477 


0.80 


1.33 


± 


0.24 


± 


0.01 


80 


± 


10 


± 


4 


1.499 


0.82 


1.23 


± 


0.24 


± 


0.01 


74 


± 


10 


± 


4 


1.522 


0.84 


0.66 


± 


0.30 


± 


0.27 


34 


± 


13 


± 


21 


1.543 


0.86 


0.57 


± 


0.29 


± 


0.11 


18 


± 


16 


± 


22 


1.565 


0.88 


0.50 


± 


0.30 


± 


0.01 


22 


± 


17 


± 


23 


1.586 


0.90 


1.18 


± 


0.35 


± 


0.01 


10 


± 


10 


± 


14 


1.607 


0.92 


1.35 


± 


0.37 


± 


0.18 


11 


± 


10 


± 


14 


1.627 


0.93 


1.11 


± 


0.37 


± 


0.14 


19 


± 


10 


± 


14 


1.648 


0.95 


1.37 


± 


0.35 


± 


0.01 


2 


± 


10 


± 


14 


1.668 


0.96 


1.82 


± 


0.43 


± 


0.01 


28 


+ 


8 


± 


12 


1.687 


0.98 


1.16 


± 


0.40 


± 


0.84 


8 


± 


14 


± 


34 


1.707 


0.99 


1.47 


± 


0.46 


± 


0.01 


11 


± 


14 


± 


21 



the Kq(1A30) resonance, and another n state. An "NR" term, defined as a constant everywhere on the Dalitz plot, 
was also included in the S'-wave 



Ci(s) 



(1430) 

Xj^(p, r R ). 



(17) 



The P- and D- waves were defined as in Eqs. (J2J and (JHJ. 

For purposes of comparison, this fit is made again, exactly as before, except that the resonance parameters indicated 
in Tableware used to replace those from Ref. p. Both the k and K£ (1430) isobars included in the 5-wavc in Eq. i|17|) 
have masses and widths that are allowed to vary. The phase convention is defined, as before, by Eq. Q. As found 
in Ref. 0, the amplitude and fraction for J4fj"(1410) are negligibly small. This resonance is, therefore, also omitted 
from this fit which is labelled the "isobar fit" . The couplings and fractions obtained are summarized in Table [n] It 
is seen that the NR term contributes modestly to the decays in this model. Its presence is, however, important as it 
interferes with the re, destructively at K~tt + threshold, not at all at 780 MeV/c 2 , and constructively at higher mass. 
Without the NR term, the S'-wave form does not fit the data well. All these results, including Breit-Wigner masses 
and widths obtained for the S- wave states, agree well, within uncertainties, with those in Ref. 
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FIG. 4: Comparison between MIPWA fit and data for D + — > K~ir + ir + decays. For each event, the distributions of (a) the 
smaller K~ir + , (b) the larger K~ir + and (c) the n + n + squared invariant masses are plotted as points with error bars. Results 
of the fit (solid histogram) are superimposed, (d) The Dalitz plot folded about the axis of symmetry resulting from the identity 
of the two 7r + mesons. The quantity plotted is the normalized residual (n/ — no)/ \jrif + n a defined in the text. The histogram 
is the distribution of normalized residual values and the curve is a Gaussian fit to this having mean —0.015 ±0.039 and standard 
deviation 0.93 ± 0.04. 



Amplitudes from this fit are plotted in Figs. OHa)-(f) where they may be compared with the MIPWA results. As 
for the MIPWA, Eqs. 10 and JHJ) are used, this time with parameters for the isobar fit in Table HTl to compute the 
magnitudes |Cl(s)| and phases 4>l(s) for the P- and D-wave, respectively, for L = and 1. Eq. (|1 Tft is used in the 
same way to compute the S-wave amplitude. Uncertainties in magnitudes and in phases are computed using the full 
error matrix from the isobar fit and values at each s, plus or minus one standard deviation, are plotted as dashed 
curves, with shading between them, in the appropriate entries in Fig. |31 

The isobar fit constrains the S'-wave magnitude and phase to assume the functional forms specified in Eq. I|17|) 
while the MIPWA allows them complete freedom. Because of the additional degrees of freedom, the latter is therefore 
able to achieve a better description of the data by a combination of shifts in the P- and £>-wave parameters, and 
in the (ck,jk) values for the S'-wave. The results presented in Figs. |2a)-(f), and in Table ILT| illustrate this. Small 
differences between the fits in parameters for if* (1680) and if.* (1430) result in relatively large shifts in the curves 
shown in Figs. GUc)-(f). These changes propagate to the S-wave. 
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FIG. 5: Moments of the K~n + helicity angle for the MIPWA fit. In each figure, events are plotted in bins of K~ir + invariant 
mass with weights equal to Pl (cos 8) as indicated. Two combinations are plotted for each K~tt + tv + candidate. Data are 
represented as points with error bars. Events are not weighted for their estimated efficiency. The MC events, treated the same 
way, are used to show the expectation for these moments from the fit, and are shown as solid histograms. 



TABLE IV: Likelihood values for fits to the E791 K~tt + system from D + -» K~ ~tv + -k + decays. The fits are described in the 
text and are labelled the same way as in Table ITU 



Model 


ln(£) 


Number of 
Variables 


NDF 


X 2 /NDF 


Probability 


MIPWA fit 


36121 


86 


277 


1.00 


47.8% 


Isobar 


36072 


16 


412 


1.08 


13.2% 


Elastic 


36092 


44 


300 


0.99 


54.9% 


Unitary 


36004 


44 


195 


2.68 


~ 



The shapes predicted by both fits for the S'-wave phase and magnitude, are shown in Figs. CH a ) and (b). Some 
differences are seen in magnitudes from K ~tt + threshold up to about 900 MeV/c 2 , and in both phase and magnitude 
above the A"q(1430) resonance. These effects are correlated with one another and with the differences in the P- and 
D-waves noted above. Similar effects arc observed in tests made on a large number of MC samples, with sizes similar 
to that of the data. Approximately 15% of these samples, generated with the distribution predicted by the isobar 
fit, give MIPWA results with similar shifts in P- and D-wave parameters, and in the associated differences in 5*-wave 
observed in the data. The MC tests are discussed in Appendix lAl 
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The significance of any differences between amplitudes obtained in the two fits is evaluated by comparing their 
abilities to describe distributions of kinematic quantities observed in the data. Plots similar to those in Figs. 01 and 
03 are made showing similar, excellent agreement between fit and data. Using the method described in the previous 
Section, the distribution observed on the Dalitz plot is compared, quantitatively with that described by the isobar 
fit results. A value for x 2 /NDF = 1.08 is obtained, and can be compared with x 2 /NDF = 1.00 for the MIPWA. These 
results are included in Table IIVI Differences between the two fits in predicted populations of bins in the Dalitz plot 
are all less than their statistical uncertainties. It is evident that both MIPWA and isobar fits are good and that no 
statistically significant distinction between these two descriptions of the data can be drawn with a sample of this size. 



VI. COMPARISON OF MIPWA WITH ELASTIC SCATTERING 



It is interesting to compare the amplitudes Cl(s) defined in Sec. IIIII and measured in Sec. IIVI with those from 
K~ir + scattering, 7l(s). The relationship between Cl and Tl is given by Eq. JfJJl. If the K~ix\ systems produced 
in D + — > /sr~7r+7r+ decays do not interact with the bachelor 7r+, then the factor Vl(s) describes the production of 
K ~ir + as a function of s from these decays. Also, under the same assumptions, the Watson theorem [24| requires 
that, in the s range where K~tt + scattering is purely elastic, 'Pl(s) for each partial wave labelled by L and by iso-spin 
I, should carry no s-dependent phase. In other words, 4>l, the phase of Cx(s) for each partial wave, should differ, at 
most, by a constant relative to that of the corresponding elastic scattering amplitude Tl(s). The magnitudes |Cx(s)| 
and |Tt(s)| could differ, however, due to any s-dependence of the production rate of K~ir + systems in D + decay. 

The validity of the Watson theorem therefore relies on the assumption that no final state scattering between (K~n^) 
and 7r+ occurs. This assumption, for decays such as those studied here in which the final state consists of strongly 
interacting particles, has often been assumed to hold. However, it has never been objectively tested. The MIPWA 
results from the present data provide, therefore, an interesting opportunity to make such a test, and also to examine 
the form for the production factor Vl{s). 



A. K 7r + Scattering 

In the S-wave, below Krj' threshold at 1.454 GeV/c 2 , K~tt + scattering in both the isospin I = h an d I — § Kir 
amplitudes is predominantly elastic. Scattering into Kr\ at a lower threshold is strongly suppressed by the 5C/(3) flavor 
coupling to this channel. This has been confirmed by the LASS collaboration in energy independent measurements 
of partial wave amplitudes for A'~7r + scattering through, and beyond this range |l9j . I = \ components of S-, P- 
and D-waves were extracted from the total using measurements of the / = | scattering from K + — > K + n + n data 
Ref. plj . Scattering in higher angular momentum waves can become inelastic at the lower Ktttt threshold. It was 
observed, however, that, in the LASS data, P-wave scattering remained elastic up to approximately 1050 McV/c 2 . 
For the D-wave, no significant clastic scattering was observed. 

In the elastic region, the I = | component of the S-w&vc K~tt + amplitude was fit, by the LASS collaboration, to 
a unitary form 

T (s) = sm[ 7 ( S )- 7o ]e^-^, (18) 
where the phase 7 = 7s + 73 + 7o is made up from three contributions: 

cot7fi = (19) 
To =0 (arbitrary offset). 

The first is a non-resonant contribution defined by a scattering length a and an effective range b. The second 
contribution r ) R has parameters tti r and r^, the mass and width of the -/^Cq(1430) resonance. In the LASS analysis, 
the arbitrary phase 70 was set to zero. The I = h P-wave and U-wavc amplitudes measured by LASS were found to 
be significant in this invariant mass range. 



B. Test of the Watson Theorem 



In Figs. Efa)-(c), direct comparisons are made, respectively, between the S-, P- and D-wave phases determined by 
the MIPWA fit to data from this experiment and the I = \ data from LASS. The 5-wave phase measurements, and 
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the curves for the P- and £>-waves resulting from the MIPWA, previously shown in Fig. EJ arc plotted, respectively, in 
Figs.|SJa), (b) and (c). The LASS measurements are superimposed, as x's with error bars, in the appropriate places 
in the figure. 




FIG. 6: Plots show comparisons between results obtained from the MIPWA, described in Sec. II VI for phases of the (a) S- wave, 
(b) P-wave and (c) -D-wave K~n + scattering amplitudes with measurements made by the LASS collaboration fTgj . Solid circles 
in (a) are the phases 7^ obtained at each invariant mass Sk from the MIPWA. The vertical dashed line at the Kr\ threshold 
(1454 MeV/c 2 ) indicates the upper limit of the range where K~ir + scattering is predominantly elastic in the S-wave. Solid 
curves in (b) and (c) enclose the zones within one standard deviation of the MIPWA P- and D-waves, computed as described 
in Sec. II VI Arrows indicate the position of the Ktvtt threshold, at which K~ir + scattering in the P- and D- waves can become 
inelastic. In (b) a further arrow at 1050 MeV/c 2 indicates the approximate invariant mass at which the P-wave data from 
Ref. fl9t were observed to become inelastic. In (d)-(f) the results from the "elastic fit" are shown in place of those from the 
MIPWA. In all plots I = \ measurements from the LASS experiment of phases for K~ir + scattering in these waves are shown 
as x's with error bars indicating statistical uncertainties. 

An obvious feature in the comparison is the overall shift in phase of the S-wave in these data relative to that in 
the LASS measurements. This feature is also evident from the examination of the K ~ir + asymmetry in the Dalitz 
plot reported in Scc.[n] Another feature of the «S*-wave comparison is that, for invariant masses near Kir threshold, 
the phases for the two sets of data show a somewhat different dependence on s. 

The P- waves also differ in the mass range from the K*(892) peak, through the region where LASS observed 
if* (1410) scattering to become inelastic, at approximately 1050 MeV/c 2 . This difference may arise, in part, from the 
paramctrization of this wave given in Eq. With more than one resonance described by Brcit-Wigncr propagators, 
this may not be unitary. The phase measured in the D-wave in this experiment agrees well with that measured by 
LASS. However, as verified by the LASS data, the scattering in this wave is no longer elastic beyond K* (892)ir + 
threshold. 

The observed shift in S'-wave phase and difference in slope, and the difference in P-wave phase behaviour evidenced 
in Figs.Eta)-(c) do not conform to the precise expectations of the Watson theorem. 



15 



C. Fit with LASS Model for S-wave Phase 



Some of the discrepancies noted above could arise from the modelling of the P-wavc. A different model could 
result in a different dependence on s of the S-wave measured here. To judge the significance of the observed discrep- 
ancies, therefore, a fit is made to the data in which the S'-wave phase is constrained to precisely follow the LASS 
paramctrization in Eqs. (|18I19() for invariant masses below Krj threshold. The mass and width of the ATq(1430) and 
the parameters a and b are required to assume the values obtained by LASS. However, the overall phase 70, all phases 
above K rf threshold, all magnitudes throughout the entire range of s, and the complex couplings for P- and D-waves 
are determined by the fit. 

This is labelled as the "elastic fit". A value 70 = (—74.4 ± 1.8 ± 1.0)° is obtained. The isobar couplings and 
resonance fractions obtained are listed in Table [D] The A"j*(1680) resonance has a more significant contribution to 
this fit than in the MIPWA. 

The phases obtained for the three partial waves from the elastic fit are compared, in Figs. EJd)-(f), with those 
measured in the LASS experiment. The comparison is shown in the same way as in Figs.|BJa)-(c) for the MIPWA fit. 
The shape of the S'-wave phase is, as required in this fit, in perfect agreement with the LASS results. The large offset 
in overall phase, 70 does, however persist. Additionally, both the P- and -D-wave phases now show larger differences 
than before. The D-wave phase shifts by ~ 50°, and the P-wave phase shows significant differences in the region 
between the if* (892) peak and the effective limit of elastic scattering at ~ 1050 MeV/c 2 . 

This fit provides another excellent description of the data, with x 2 /NDF = 0.99 and probability 55%, as recorded in 
Table Hvl This is comparable with both the isobar and MIPWA fits. However, if these observations are predominantly 
of production of / = | K~ir + systems, the phase variation required by the Watson theorem, is not observed in these 
data. 



D. Production Rate for the K ix + System 

For purely elastic scattering, the Tj, amplitudes are required to be unitary, as given by Eq. <|18|l . Introducing this 
into Eq. JBJ leads to 



(vdsy 

sin[ 7 (s) -7o]e l[7(s) - 701 . (20) 



cm - ( f j {^i ) * 



For L = 



| Po(8) | . I JL) , (21) 



sfs) Vsin[7(s) -7o] 

Structure in the s-dependence of the S-wave magnitude, Co(s) can thus come either from the phase 7(s), from Vq{s), 
or from both. It is of interest to study these possibilities, and to see if the data can be described by a unitary 
amplitude, in which Pq(s) would be independent of s. 

The data from the MIPWA are examined to see if the S-wave can be described by a unitary amplitude, such as 
that given in Eq. (|20() . Setting L = and Vq(s) = P (a constant), a value for 70 is determined by minimizing the 
quantity 



\Vo(s k )\-P 



where |"Po(sfc)| are computed from Eq. (|21|l . for the values of S-wave amplitude, Co{sk) = cj.e' 7fc , determined by the 
MIPWA fit, and cr(Vo) are the associated uncertainties. The summation in Eq. I|22|) is made only for the Nk v ' values 
of s k up to the K~r]' threshold. The value 7 = (-123.3 ± 3.9)° is obtained, with P = 0.74 ± 0.01 (GeV/c 2 )~ 2 . 
Fig. 01a) shows that this value for 70 is approximately equal to the measured S-wave phase at K~ir + threshold, 
consistent with the physical meaning of this parameter in the formulation in Eq. (|20|l . 

Inserting this value for 70 into Eq. |JHJ, the quantities l^oSfcl are plotted in Fig. [7] The solid, horizontal line indicates 
the value for P obtained from the fit. The points are seen to lie close to this line, showing very little dependence 
on s in the invariant mass range from K~tt + threshold up to about 1.25 GeV/c 2 . From 1.25 to 1.5 GeV/c 2 , strong 
variation is observed. In this region, as seen in Fig.|S{a), the value of sin(7 — 70), which appears in the denominator 
of Eq. (|21f) . is approximately zero. 
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FIG. 7: The quantities p/-\/(sfc) x |Co(sfc)| x JF^,/ sin(7k — 70) plotted as solid circles for each point obtained for the .S'-wave 
amplitude in the MIPWA fit described in Sec. II VI Three points between 1400 and 1450 MeV/c 2 are omitted from the plot as 
their values for sin(7fc — 70) are very small. Their values are either off-scale or their errors extremely large. The region between 
the dashed lines shows the one standard deviation limits of this quantity for the S'-wave amplitude obtained from the isobar 
fit. The inset shows, as small open circles, the quantities |Tb(s)|/ sin(7 B + 7^) taken from the LASS experiment. The Krf 
threshold is indicated by dashed, vertical lines in both plots. 



Also shown in Fig.[7|is | T^o ( s ) I for the isobar fit. The region between dashed lines corresponds to the one standard 
deviation limits for this quantity, computed from Eq. (|21[) with the same value of 70 as used above. Values for 
magnitude and phase of the S'-wave amplitude, and their statistical uncertainties, are computed from Eq. I|17f) with 
parameters and error matrix from this fit. The behaviour of |"Po( s )l derived from the isobar model fit matches that 
observed in the MIPWA points well. 

The inset in Fig. [7] shows the corresponding quantities |To(s)|/ sin(7 B + 7 R ) for the points measured for K~ir + 
scattering in the LASS experiment. From Eqs. (|18|) and (|19|) it is seen that this is expected, in the range up to K~r)' 
threshold, to be unity. It is seen that this condition is met by the LASS data. 

It is concluded that the factor IT-'o (s) | in Eq. (|21fl that describes production (and possible re-scattering) for K~ir + 
systems in the D + decays examined here, shows little dependence on s up to about 1.25 GeV/c 2 . At this point, a 
significant dependence on s is seen. This behaviour is qualitatively different from elastic scattering. 



VII. SYSTEMATIC UNCERTAINTIES 



The major source of systematic uncertainty in the MIPWA results arises from the difficulty, with a sample of this 
size, of reliably characterizing the structures, other than the if* (892) resonance, in the reference waves. To estimate 
this effect, a large number of samples of MC events, each of which is of a size similar to the data (~ lhK events) 
reported here, are examined. These are generated with the parameters determined by the isobar model fit described 
in Sec. |3 with the backgrounds best matched to the E791 data. Each sample is subjected to a MIPWA fit, and the 
differences between generated and fitted values for S-wave magnitude and phase at each of the 40 invariant masses 
are examined. For most samples, fits obtained match the isobar model well. Variations in the significance of the 
if* (1680) and K^iUSO) sometimes lead to variations in the reference waves that propagate to distortions in the 
S-wave solutions found. These tests provide estimates of systematic uncertainties for the S'-wave magnitudes that 
range from ~ 50% of the statistical uncertainty, for invariant masses below 800 McV/c 2 , to an insignificant level 
for higher masses. For the S-wave phases, the systematic uncertainties are found to average ~ 72% of the statistical 
uncertainty. 
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The second largest uncertainty arises from the smearing of events near the high mass boundary of the Dalitz plot 
which results from the resolution in three-body mass M. This directly affects part of the if j* (1680) band. Events in the 
region of M closest to the D + mass are fitted separately, and the results compared with that from the larger sample. 
Average systematic uncertainties arising from the effects of smearing are determined to be 7% of the statistical 
uncertainties for magnitudes and 14% of the statistical uncertainties for phases. Other effects are studied. These 
include the uncertainty in precise knowledge of the background level, variations in the values assumed for the radii 
r R and r n , or for the mass and width for the if* (1680) resonance. All these other effects are found to be small. 

A further source of systematic uncertainty arises from variations in the presumed resonant composition of the 
K~ir + P- and £>-wavcs. The if* (892) resonance obviously contributes, and it is clear that a contribution from 
a higher resonance, must also exist. What is less clear is the identity of this resonance - if*(1430), 7^(1680) or 
for if*(1410). Fits are made with various combinations of these resonances. It is found that systematic shifts are 
negligibly small in most cases. Fits where only if* (892) and if.* (1430) are included do lead to shifts comparable to 
the statistical uncertainties in the lowest five magnitudes. At higher invariant masses, the effects become smaller. 
The phases are almost unchanged, however. 

These uncertainties are combined in quadrature and listed for each invariant mass in Table ITTll and for the reference 
wave parameters in the MIPWA fit in Table ITT1 

VIII. SUMMARY AND CONCLUSIONS 

A Model-Independent Partial Wave Analysis (MIPWA) of the S-wave K~tt + system is made using the three body 
decay D + — > K~it + it + . This is the first time such a technique has been used in studying heavy quark meson decays, 
and new information on the K~ir + system is obtained, including the invariant mass range below 825 McV/c 2 . The 
isospin I of the S- wave measured is unknown, and the P- and D-w&ves are assumed to be / = h. It is possible to 
modify these assumptions, provided independent information on the / = | components is available. The method does 
not assume any form for the energy dependence of the S-wavc. However, it does so for the P- and D- reference waves. 
The P-wave is described as the sum of a Breit-Wigner propagator term for the if* (892) resonance, and a similar 
term, with a complex coefficient, for the if* (1680). The if*(1410) is found to have an insignificant contribution to 
the decays, and is omitted from this wave. The D-w&ve is described by a single Breit-Wigner term for the if | (1430) 
resonance, with a further complex coefficient. The results obtained in Fig. Eland Table ITTT1 depend on the accuracy of 
this description of the reference waves. 

Results of the MIPWA are compared with a description of the S-wave amplitude that includes Breit-Wigner 
K, ifp(1430) isobars and a constant, non-resonant (NR) term similar to the description used in Ref. [||. At the 
statistical level of this experiment, differences between the MIPWA and the isobar- model result are not found to be 
significant, and both provide good descriptions of the data. A closer examination of the phase behavior in the low 
mass region below 825 MeV/c 2 , the limit of measurements of if~7r + elastic scattering from the LASS experiment 
[lflj . is of considerable importance to the further understanding of scalar spectroscopy. The data here provide new 
information in this region, but the error bars are large compared to those typical for the LASS data. We note 
that, since these data became available, a fit that includes requirements of chiral perturbation theory has been made 
together with the LASS I = 1/2 measurements and data from J/ip — > if *(892)if ~7r + . This fit finds a k pole at 
(740lgg) — i(342 ± 60) MeV/c 2 j34j. A full understanding of scalar K* spectroscopy may, nevertheless, need to wait 
until larger data samples become available. A better consensus on the proper theoretical description of such states 
and the need for, and the form of, any accompanying background amplitudes may also be required. 

The phases observed in the S- and P-waves do not appear to match those seen in the I = i elastic scattering in 
reference 0] . The P-wave phase does agree well. Constraining the energy dependence of the S'-wave phase to follow 
that observed in (I = |) 2f7T elastic scattering, in the range where s lies below if 7/ threshold, does lead to a good fit 
to the data. However, an overall shift in phase of (— 74.4 ± 1.8 ± 1.0)° relative to the P-wave is still required. This 
constraint also results in a shift of approximately —50° in the D-wave phase. It also makes agreement in P-wave 
phase worse. These results do not conform, exactly, to the expectations of the Watson theorem which would require 
phases in each wave to match, apart from an overall shift, those for if ~n + scattering for invariant masses below if 77' 
threshold. The theorem is expected to apply in kinematic regions where secondary scattering of the Kir system from 
the bachelor pion can be neglected. It is possible that, in this case, such scattering cannot be neglected, or that the 
if _ 7r + systems in D + decay are not predominantly (i = |) [35j . 

It is also found that, with a choice of phase at if ~ir + threshold relative to the P-wave 70 = (—123.3 ± 3.9)°, quite 
consistent with that measured in the MIPWA, if _ 7r + systems produced from D + decays are described well by a 
unitary amplitude (with constant production) up to a mass of about 1.25 GeV/c 2 . In this region, therefore, structure 
observed in the S'-wave magnitude is mainly associated with the variation in phase with respect to s. the invariant 
mass squared in the if~7r + system. Above 1.25 GeV/c 2 , the production rate grows, depending significantly on s. 
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The reason for this behavior is unknown. The growth observed at 1.25 GeV/c 2 could result from a significant / = | 
contribution or from re-scattering of the produced K ~ir + system and the bachelor ir + . 

The MIPWA analysis of the three-body decay of a heavy quark system described here has three main limitations. 
The first results from the way the reference P-wave is described in Eq. J7J. Using Breit-Wigner resonance forms for 
more than one resonance in the wave can lead to problems in the regions where the resonance tails dominate. The 
second limitation comes from the ability to resolve the structure in the if (1680) region properly at the statistical 
level of the E791 data. This problem may be specific to the channel discussed here, and to the particular data sample 
used. The third limitation is the lack of knowledge on any I = | components in the system. 

The first two limitations should be mitigated when much larger data samples are available. A better formulation 
for the P-wave could be to use a if -matrix, requiring more parameters. Alternatively, the P-wave, too, could 
be parametrized like the S'-wave, in a model-independent way. The third limitation can be improved when larger 
samples of K + tt + or K~ir~ systems can be studied to better understand these waves. 

Systematic studies of various heavy quark meson decays in future experiments (BaBar, Belle, CLEO-c, and hadron 
colliders), with much larger samples, may be able to use a similar MIPWA technique, with some of these improvements, 
to shed further light on important questions in light quark spectroscopy, the realm of applicability of the Watson 
theorem, etc. For studies that require an empirically good description of the complex amplitude in three bod y de cays, 
for example in the extraction of the 7 CP violation parameter recently reported by BaBar and Belle 0, l9lTi"o| , this 
technique may also be particularly useful. 

In the mean time, theoretical models of the S-wave amplitude can be compared to the data of Table III. 
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APPENDIX A: LIMITATIONS AND TECHNICALITIES OF THE METHOD 
1. Quality of S- wave Measurements 

The MIPWA analysis of the S-wave component in the observed D + — > if~7rj~7r+ decays relies upon a good 
description of the reference P- and P-waves. 

The P-wave defined in Eq. J7J) is a combination of two Breit-Wigner's (the K ^(1410) is neglected), with a complex 
coupling coefficient (two parameters). The peak regions, 0.50 < s < 0.9 (GeV/c 2 ) 2 and s > 2.2 (GeV/c 2 ) 2 , are well 
described by this parametrization, since data in these regions are likely to be dominated by these resonances. In the 
tail regions, the Breit-Wigner may be a less appropriate description of the data, since other P-wave contributions, 
from non-resonant (or I = |) sources, for example, could become more significant. 

Two regions in the P-wave where the tails of the BW's dominate are s < 0.50 (GeV/c 2 ) 2 and 0.9 < s < 
2.2 (GeV/c 2 ) 2 . In each of these ranges the P-wave is constructed from a linear combination of two small, com- 
plex numbers, one from each of the two BW tails. Both the phase and magnitude of the resultant are particularly 
sensitive to variations in the complex coupling parameters, and may not represent the P-wave well. 

The P-wave is defined in Eq. (JSJ), for this analysis, as a single, L = 2 BW function. Non-resonant contributions 
are not expected to be significant, and interference from tails of a second resonance are absent. Eq. ©, therefore, 
provides a relatively good description of the P-wave. 

Eqs. J2)l and JJJl show that the amplitude for the decays examined in this paper is a sum of six terms. Let these be 
labelled S A , V A and T> A (the S-, P- and P-wave, respectively, in the s A channel) and S B , V B and T> B (these waves 
in the s B channel). The MIPWA process extracts magnitude and phase information about the .S-wave S A from its 
observed interference with the complex sum of the other five amplitudes. 

T A =T A + V A +S B +V B +V B . (Al) 

The results expected from measurement of S A can, therefore, be characterized by the dominant terms in T A with 
which it interferes, and these depend on location on the Dalitz plot in Fig. ^ 
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FIG. 8: Comparison between fitted S'-wave amplitudes (solid circles with error bars representing statistical uncertainties) and 
generated amplitudes, represented as shaded regions between dashed curves computed as described in Sec.|3 Generated curves 
follow the isobar fit model described in Sec.0 The figure shows results from the first six (of 100) samples used in the test. 



As an illustration, consider measurement of S A in the range 1.1 < s < 2.9 (GeV/c 2 ) 2 . Here, T A is dominated by 
the if* (892) resonance band in V B (the cross-channel). Good measurements are, therefore, expected in this range. 

Next, consider the if* (892) peak region 1.5 < s < 0.9 (GeV/c 2 ) 2 . This region has T A dominated by the if* (892) 
peak in V A (the direct-channel). So good measurements are expected here too. A similar conclusion can be drawn 
for the ^(1680) peak region s > 2.9 (GeV/c 2 ) 2 . 

Relatively poor measurements are expected for the other regions since, in these, T A is not dominated by any one 
source, and is defined by a linear combination of several Brcit-Wigncr tails. So T A in these regions has phase and 
magnitude that are sensitive to the complex couplings of if J 1 (1680) and if.* (1430) resonances. 

These observations are supported by the results of the MIPWA fit shown in Fig. E{a) and (b). It is seen that 
uncertainties are small for 1.1 < s < 2.9 (GcV/c 2 ) 2 and large for s < 0.9 (GeV/c 2 ) 2 , improving towards the high end. 
In the intermediate region, 0.9 < s < 1.1 (GeV/c 2 ) 2 , the 5-wave magnitudes and phases determined in the fit exhibit 
significant deviations from the general trends of the neighboring points. 



2. MC Studies with Isobar Fit 



The if * (1680) and if* (1430) resonances represent small contributions to the Dalitz plot, and statistical uncertainties 
in their complex couplings are large enough to affect the P-wave phase, especially in the region between if * (892) and 
1^(1680), as discussed in Sec. EJJ This can lead to systematic uncertainties in the S'-wave amplitudes measured. 
MC studies are required to estimate such effects. 

MC samples of the approximate size of the data presented in this paper are generated as described by the PDF given 
in Eq. 1|14|) . Parameters from Table ITU for the isobar fit described in Sec. are used for this purpose. Background 
events whose distributions are given in Eq. i|13|) are also generated to match those thought to be present in the data. 
Events are selected according to the efficiency e(s A , s B ) across the Dalitz plot. 

Each sample is subjected to the MIPWA fit described in Sec. IIVI In Fig. S-wave amplitudes determined in 
the MIPWA for the first six of the 100 samples studied are compared with those used to generate the events. The 
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amplitudes generated come from the isobar fit, and are shown, as usual, as shaded regions between dashed curves. 
Phases are shown on the right and magnitudes on the left. Plots similar to Fig. [3] appear often. 

The S-wave amplitudes (cfce I7fc ) obtained are compared with those generated and, for each k the normalized residuals 
are used to determine systematic uncertainties discussed in Sec. IVIII 



3. Other Solutions 




FIG. 9: A second solution for the S-wave amplitude from MIPWA fits to D + — > K~ir + ir + decays with P- and D-wave 
parametrized by the k model described in the text. Plots show the (a) phase and (b) magnitude for solution B for the S- wave 
obtained by using different starting values for the amplitudes. The dashed curves delineate the regions that lie within one 
standard deviation of the isobar model fit described in Sec. El The P-wave is shown in (c) and (d) and the D-wave in (e) and 
(f). 

The fitting procedure allows a great deal of freedom to the S-wave amplitude. Consequently, ambiguities in solutions 
are anticipated. To study possible ambiguities in the MIPWA solution, fits with random starting values for the Cj,7j 
parameters, and also with different K~tt + mass slices are made. One other local maximum in the likelihood is found, 
and this is labelled solution B. The solution described in Sec. IIVI and shown in Figs.|3{a) through (f), is labelled, for 
contrast, solution A. Solution A is the only one with an acceptable x 2 /NDF and has the greatest likelihood value. So 
it is emphasized that solutions A is, in fact, unique. 

Solution B is shown in Fig. [HI It provides a qualitatively reasonable description of the distribution of the data on the 
Dalitz plot. However, this solution clearly exhibits retrograde motion around the unitarity circle as K~n + invariant 
mass increases. This violates the Wigner causality principle j^B], thus eliminating it from further consideration. 

The possible existence of other maxima in the likelihood, when all 5-wave magnitudes and phases are free param- 
eters, cannot be completely ruled out. However, the solution in Sec. IIVI is unique in that it is the only one giving an 
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acceptable fit probability. 
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